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Abstract. As an undergraduate exercise, in an article [Am. J. Phys. 80 (2012), 
708-14], quantum and classical uncertainties for dimensionless variables of position 
and momentum were evaluated in three potentials: infinite well, bouncing ball, and 
harmonic oscillator. While original quantum uncertainty products depend on h and 
the number of states (n), dimensionless approach makes comparison between quantum 
uncertainty and classical dispersion possible by excluding h. But the question is 
whether the uncertainty still remains dependent on quantum number n. In the above 
mentioned article, there lies this contrast; on the one hand, dimensionless quantum 
uncertainty of potential box approaches classical dispersion only in the limit of large 
quantum numbers (n —^ oo) - consistent with the correspondence principle. On the 
other hand, similar evaluations for bouncing ball and harmonic oscillator potentials 
are equal with their classical counterparts independent of n. This equality may hide 
the quantum feature of low energy levels. In the current study, we’ve changed the 
potential intervals in order make them symmetric for the linear potential and non- 
symmetric for the quadratic potential. As a result, it is shown in this paper that the 
dimensionless quantum uncertainty of these potentials in the new potential intervals 
is expressed in terms of quantum number n. In other words, the uncertainty requires 
the correspondence principle in order to approach the classical limit. Therefore, it 
can be concluded that the dimensionless analysis as a useful pedagogical method 
does not take away the quantum feature of n-dependence of quantum uncertainty 
in general. Moreover, our numerical calculations include the higher powers of position 
for potential. 
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1. Introduction 

Since the advent of quantum mechanics, there have been substantial works attempting 
to approximate the novel quantum notions to the old classical concepts. Among them, 
we can refer to correspondence principle [1], bold classical path [2], quantum Virial 
theorem [3], joint distribution function of position and momentum [4], Dirac’s analogy 
[5], etc. Such theoretical researches are still referred to by scholars as major sources on 
the foundations of quantum mechanics, among which we can mention studies [6, 7, 8]. 

On the other hand, some studies have attempted to build a probabilistic structure 
for inherently predetermined variables of classical mechanics in order to compare 
essentially statistical quantum observables. According to Curtis and Ellis [9], an 
introductory physics course utilizes a historical Newtonian approach which is identihed 
by instantaneous values for position, speed and acceleration, whereas a quantum course 
is featured by conceptually probabilistic observables. 

“Although the traditional tracking of instantaneous positions (like a series of 
snapshots) has conceptual advantages, the use of position probabilities (like 
a time exposure that reveals motion through the degree of overexposure) has 
other advantages. It allows macroscopic and microscopic objects to be studied 
by similar techniques. It is easily extended to include many-body interactions, 
since the probability distributions can be superimposed. ... Another advantage 
of the probabilistic approach is that it provides a convenient way of performing 
numerical computations for potentials that do not have an analytic solution. ” 

This approach leads to the comparison of classical and quantum probabilities in high 
energies by applying the correspondence principle. In other words, quantum distribution 
wiggles so rapidly in a scale set by the classical amplitude that only its mean can be 
detected at these scales, and this agrees with classical probability [10]. Hence, squared 
wave function and classical probability density, if averaged in a hnite interval, seem to be 
indistinguishable for large quantum number n of one-dimensional problems. In a study 
conducted by Robinett [11], this point of view was demonstrated by several examples. 

Here we highlight the structure of classical probability in more details. Basically, 
position probability density of one point particle (or an ensemble of identical particles) 
with its motion equation in one dimension x = X{t) can be written as 

p{x,t) (X 6{x — X{t)). (1) 

On the other hand, for a time-independent potential, we have p{x,t) = p{x)T{t) and 
therefore, we can track time-independent position probability density p{x) via the 
probability of Ending the particle in a small region, dx, proportional to the amount 
of time, dt, the particle elapse there. 

probability[(a:, x + dx)] = p{x)dx, 

p(x)dx X dt ^ nix) x , , oc — , =, 

vix) ^E-Uix) 


( 2 ) 
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, where v is speed in terms of x, and E and U stand for total and potential energy, 
respectively. The equivalency of equations (1) and (2) concerning position can be 
reached by the property of delta function, as can be seen in the following equation; 

Six-X{t)) = , ti = X-\x). (3) 

For time-independent potentials, we have a single-valued function for speed v{x) in 
terms of x. Therefore 


vix) Six - Xit)) = '^Sit-ti) , vix) = \X'iti)\. (4) 

i 


Now, substituting (1) on the left side of (4) and integrating in respect to dt on both 
sides of (4) result in (2). As an example, for a Harmonic oscillator A’(t) = Asin(a;t±(/)) 
and vix) = Acucos(sin“^-|) = A? — hence R = ;^[sin“^ is unique so that we 

do not have summation on the right side of (4), unlike the free fall equation of motion, 
which leads to two terms for this summation. 

We can also derive similar equations for the classical probability density of 
momentum. For instance, an analogous equation for (2) in the momentum space would 
be [11] 

pip) oc (5) 


, where F refers to the force. 

Equation (4) can be interpreted as a classical formula where determined position 

and speed for point particle are matched together. However, knowing this determination 

can not prevent us from calculating classical dispersion as (Ap)^j by means of 

averaging in respect to classical probability densities, so that it can be compared with 

quantum uncertainty relation (Ap)^^^. 

As Hamiltonian approach can substituted for Newtonian mechanics, it can also be 

applied for the classical probability of the phase space in case we write a delta function 

2 

on the subject of Hamiltonian H = ^ + U (x) in a given Energy 


pix, p) oc (5 


V 

2m 


+ Uix) - E 


( 6 ) 


, where constant F in a given potential, directly refers to one particle or an ensemble of 
identical particles. 

Straightforwardly, classical averaging can be obtained just in terms of position 
probability density [12] (see Appendix A) 


{fix, P)) = o / da;p(x) [/ (x,-^j2m[E -Uix)]) + f (x, ^j2m[E -Uix)] 


( 7 ) 


Devi and Karthik [12] have chosen dimensionless variables of position and momentum 
as 


X 


X / X^uifix , 


P P/Pmax) 


( 8 ) 
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, where Xmax and Pmax are classical return points; then they analytically compared 
classical dispersion (AX)^j (AP)^j with quantum uncertainty of (AX) (^-P) for 
three potentials of infinite well, bouncing ball, and harmonic oscillator. 

The following dimensionless classical position and momentum dispersions 

(AA-)^, = (A^>^, - (A)^„ 

= (^’^> 0 , - (P)l 

could be evaluated using (7), whereas dimensionless quantum uncertainty could be 
evaluated by averaging on squared wave functions, so Xmax{n) and Pmax(ri) are derived 
as functions of n from Eqm = Pd- 

With some generalization, their results for three mentioned potentials are as follows: 


• Symmetric or non-symmetric potential box, where A is a constant 

f 0 <x<Aoit)<x<A 

1 oo otherwise 

Imi (aa)' (ap)' = (AA)= (AP),^, = | 

n^oo \ /qmV /qm 1/12 

• Bouncing ball| (non-symmetric linear) potential 

{ oo for X < 0 

, 

ax for X > 0 

(aa)" (ap)' = (AA);^, (AP)^, = 

• Harmonic oscillator (symmetric quadratic) potential 
U (x) = ax^, 

(aa)' (ap)' = (AA);qAP)^, = 

, where a is a positive integer constant which is neutral when choosing dimensionless 
variables (just like m and h in the relevant Schrodinger equation of Appendix B). 

Evaluations of potential box (10) approach classical dispersion only in the limit of 
large n, which is consistent with the correspondence principle. Namely, its dimensionless 
quantum uncertainty depends on quantum number n. But for bouncing ball (11) and 
harmonic oscillator (12) potentials, quantum uncertainties are independent of n and 
exactly equal to classical evaluations. The reason why the quantum uncertainty of 
bounding ball and harmonic oscillator potentials are to be independent of n whereas 
that for the potential box is not (or further, for other powers of x in the potential 
except for powers 1 and 2) could be left as an exercise, [hint: use < x >oc< 7/ > or 
< x^ >oc< U > OT < lA >oc< 7/^ > and consider Virial theorem]. 


symmetric 

non-symmetric 


( 10 ) 


135 


( 11 ) 


1 

4 


( 12 ) 


f The term “bouncing ball potential” describes a uniform gravitational field above a rigid flat surface 
if ax —>■ mgz in (11). Quantum solutions for bouncing ball potential are just the odd states of the 
symmetric linear potential. For more information about this potential, see [11, 13, 14]. 
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We could obtain uncertainty relations for original variables of x and p instead of 
dimensionless equations (10), (11), and (12). Original uncertainty products 

(Ap)^„ > f (13) 

depend on h] also it depends on n since we have integrands as 'ipn{x) and we take 
integrals in respect to variable x] then n is preserved. For example, we can show that 
(A£)qj^ (Ap)q^ = (n + l/2)^h^ for the symmetric harmonic oscillator [15] [hint: use 
properties of Hermite functions or equivalently use the effect of annihilation and creation 
operators on a given state]. When we chose dimensionless variables of (8) at the left 
side of (13), it causes the dimensional constant of h (or h) to be prevented at the right. 
Although excluding h is an advantage of the dimensionless approach of uncertainty 
and makes its comparison with classical dispersion possible, regarding n-dependence as 
an other quantum feature of quantum uncertainty, there is still the question whether 
this dimensionless uncertainty remains dependent on quantum number n under other 
circumstances. 

In the present study, we bring into consideration dimensionless quantum uncertainty 
versus dimensionless classical dispersion, as obtained by averaging on classical 
probability density. First of all, we numerically calculate dimensionless classical and 
quantum uncertainty for symmetric and non-symmetric potentials from powers of 
coordinate. Then, we narrowly focus on linear and quadratic potential solutions and 
compare them with the results of the [12] (equations (11) and (12)). The key question is 
whether the dependence on n for the dimensionless uncertainty of these two potentials 
can be retrieved if their potential intervals be changed. To seek an answer to this 
question in section 3 of this article, we will focus on the uncertainty of symmetric linear 
potential and compare it with that of the non-symmetric linear potential (bouncing ball) 
of (11). Also, in section 4, we deal with the uncertainty of the non-symmetric quadratic 
potential and compare it with that of the symmetric quadratic potential (harmonic 
oscillator) of (12). In the above mentioned sections, regarding linear and quadratic 
potentials, it is shown that although there is a potential interval in which dimensionless 
quantum uncertainty is equal with its classical counterpart -independent of quantum 
number n, according to [12], another interval can be chosen to reveal the dependence on 
n. Therefore, this change can be considered as an emphasis on dimensionless approach 
as a practical pedagogical method that does not take away the quantum feature of 
n-dependence in general. Moreover, the solutions of potentials with higher powers of 
coordinate conhrm n-dependence behavior in dimensionless uncertainty. Last but not 
least, we discuss shortly the quantum uncertainty concept being preserved in the classical 
limits and also, look at the effect of measurement on the interpretation of the classical 
limit of a non-degenerate wave function as a classical ensemble. This study aims to 
emphasize on both numerical and analytical curricula. Accordingly, the first and third 
appendices contain analytic contexts used in the main text and the second appendix 
includes some basic materials used for numerical solutions of Schrodinger equation in 
the matrix representation method. 
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(c) (d) 

Figure 1. Classical (dashed) and quantum (solid) uncertainty of dimensionless 
variables for symmetric potentials of equation (14) (horizontal axis stands for powers 
of |a;|). 1(a) n=0, 1(b) n=l, 1(c) n=2, 1(d) n=10. 


We consider symmetric potentials 

U{x) = a\x\'^, where 6 = 1, 2,..., 10 (14) 


and non-symmetric potentials 


U{x) = 



for x < 0 , 
for x > 0 


(15) 


, respectively, in the Schrodinger equation, 

r|^ 7 /; ( '7*') 

- + U{x)i:{x) = Ef^ix) (16) 

and numerically calculate the amount of quantum uncertainty and classical dispersion of 
each potential for dimensionless position and momentum variables dehned in (8). These 
dimensionless quantities are bare numbers such that |X|, \P\ < 1 in a bound potential 
[ 12 ]. 
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0.09 I-1-1-1-1-1-1-1-r 


0.09 I-1-1-1-1-1-1-1-r 



(c) (d) 

Figure 2. Classical (dashed) and quantum (solid) uncertainty of dimensionless 
variables for non-symmetric potentials of equation (15) (horizontal axis stands for 
powers of x). 2(a) n=0, 2(b) n=l, 2(c) n=2, 2(d) n=10. 


At first, eigenvalues {Eqm{n)) and eigenfunctions ('0„(a:)) of each Schrodinger 
equation in a given potential are derived and then the related ( AX j ( AP j is 

evaluated by averaging on each state in the given potential. Next, the classical dispersion 
(AA)^j(AP)^j is calculated using (7) with Ed = Pqm, and compare it with each quantum 
uncertainty. 

Eigenvalues and eigenfunctions can be derived easily by creating a matrix 
representation for the second derivative operator at hrst (see Appendix B); then the 
related arrays of potential (as a diagonal matrix) are added to the arrays of the second 
derivative matrix and hnally, the eigenvalues and eigenvectors of the summation matrix 
can be taken using a simple assigned command for this purpose in computer programs. 
Some useful Matlab codes for educational quantum problems are created in [16]. 

Quantum uncertainties and classical dispersions for symmetric potentials of (14) 
are shown in hgure 1 in different quantum states. Also, hgure 2 shows the results of 
non-symmetric potentials of (15). 
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Since b = 2 in (14) refers to the symmetric harmonic oscillator potential, quantum 
uncertainty for b = 2 in hgure 1 is equal to the analytic result of (12), or it is exactly 
equal to classical dispersion for this potential -independent of the quantum number n. 
Also 6 = 1 in (15) refers to bouncing ball potential and again, its quantum uncertainty 
(see & = 1 in the hgure 2) is exactly equal to its classical counterpart, which is in 
agreement with (11) -independent of n. 

In the next two sections, we will argue about 6 = 1 in hgure 1, which belongs to 
numerical evaluations for symmetric linear potential 


U{x) = a\x\ — A < X < A 


( 17 ) 


with amplitude A, and 6 = 2 in hgure 2 is related to numerical evaluations for non- 
symmetric harmonic oscillator 


U(x) 


oo for a; < 0 , 

ax'^ for a; > 0 . 


(18) 


3. Comparison between quantum and classical uncertainties for symmetric 
linear potential 


Bouncing ball potential is non-symmetric; however, by considering the symmetric 
potential of U = a|a:|, we can observe that at least for n = 0 (5 = 1 in hgure 1(a)), 
classical and quantum uncertainties have diherent values. In fact, bouncing ball 
potential results in half of the wave functions of symmetric linear potential; in the other 
words, we will have solutions with "^(0) = 0 (odd wave functions). On the other hand, 
symmetric linear potential inserts all of "^(0) = 0 and "^'(0) = 0 solutions § and quantum 
uncertainties in diherent even n become diherent from each other quantitatively. Despite 
this, according to the correspondence principle, in the large n, uncertainties become 
approximately equal for both odd and even solutions. 

Considering quantum averages, we seek for averages which cause diherence between 
classical and quantum uncertainties of the linear potential. We know that is 

Also, it can be derived that and are 

remains zero since its integral in the position 

in bounds and for the second moment {P'^^i 

can derive (p^) = aE{n) and {P'^'^ oc (p^)/P from the quantum Virial theorem; then 

^P^^ is independent of n). Therefore, the diherence between classical and quantum 
uncertainties for 6 = 1 in hgure 1(a) must be created by n-dependence of the second 
moment of position Figure 3 displays quantum and classical (W^) for the 

symmetric linear potential in terms of n and explicitly shows that ~ 0.72 

\ / n=0 


zero for the symmetric potentials, 
independent of n (the hrst moment 
space depends on the value of fj^x] 


§ Matrix elements of the quantum bouncer [14] do not involve even wave functions of symmetric linear 
potential or those for which = 0. 
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(a) (b) 

Figure 3. Classical (dashed) and quantum (solid) belonged to solutions of 

symmetric linear potential in terms of energy levels n (X is a dimensionless position 
variable defined in Eq. (8)). Figure 3(a) shows clearly the value for n = 0, while in 
hgure 3(b) the same graph is vertically magnified to show other states except n = 0. 


keeps distance from that of other states. We have magnified vertical axis scale for 
better visibility of other states in figure 3(b), such that 0.72 could not be seen. Difference 
between quantum uncertainty and classical dispersion at the lowest quantum state of 
symmetric linear potential (6 = 1 in hgure 1(a)) arises from the difference between 
0.72 in hgure 3(a) with 8/15, which is related to the classical amount and the amount 
calculated for wave functions of bouncing ball potential [14] (odd wave functions of 
linear potential). It can be observed that the diherence with classical limit for other 
even states in hgure 3 is slight and shows decrease. So, this diherence is not dramatic 
in hgure 1 for the next even n. 


xn >(xn >(xn ...~(xn =(xn ={x^' 

n=0 \ / n=2 \ / n=4 \ / n=l \ / n=3 \ / n=5 


and 

' = /xA = (x2) 

n=2k \ / n=2k-\-l 

Therefore, uncertainty becomes dependent on n. 


lim (X‘ 

n—>-OD 


cl 


k = 0,1,2, 


(19) 


( 20 ) 


4. Comparison between quantum and classical uncertainties for 
non-symmetric harmonic oscillator 


Wave functions are for 


U{x) = 


( 21 ) 


oo for X < 0, 

ax^ for X > 0, 

odd solutions of the harmonic oscillator potential located in the interval [0, oo) with 
doubled normalization constant. So 



non—symmetric 



^ symmetric 


1 

2 ’ 
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P ) =( P 

non—symmetric 


symmetric 


p^ 


= P 


non—symmetric 


/ symmetric 


= 0 ; 

1 

2 ‘ 


( 22 ) 




Non-symmetric can be evaluated analytically, as shown in Appendix C. In the 

absence of even wave functions, odd wave functions here are starting from n = 0; 
therefore, 2n + 1 in (C.15) is transformed to n here. 

2 (2n + l) {2n)\ 


x\ 


/ (2n+l)^-n 71-2 + 3)2 2^” (n!) 

As a distinctive case, for n=0, we have 

2 1 


2 • 


a: 


n=0 




which results in 


aa:) 

/ qm 


AP) 

/ qm 


n=0 


1 4 \ 1 

-- ~ 0.0378 . 

2 37r/ 2 


(23) 


(24) 


(25) 


This is in agreement with fe = 2 in hgure 2(a). Thus, uncertainty is again dependent on 


n. 


5. Discussion 

In this section, we briefly take a look at two general points which should be considered 
when studying classical limits. 

5.1. Quantum uncertainty in classical limit 

In addition to topics of sections 3 and 4 as the main goals of this article, hgures 1 and 
2 also show that position-momentum quantum uncertainty relation does not vanish in 
classical limit. In other words, quantum uncertainty tends to classical evaluations in 
the classical limit, but not to zero. Namely 

lim (ax)' (ap)' =(AX)l(AP)lTt0 

n—>-oo V / qm V / qm 

Therefore, in the classical limit of n —?■ oo, quantum particles still have quantum 
uncertainty in its quantum context. It means that the equality relation between n ^ oo 
limit of quantum uncertainty and classical dispersion is approaching averages, not 
approaching concept of quantum uncertainty to that of the classical. 

5.2. Ensemble’s interpretation in classical limit 

There are articles such as [6, 17] that emphasize that the classical limit of a wave 
function describes an ensemble of classical particles -not an individual particle. For 
example, Ballentine et ah [6], when challenging Ehrenfest theorem, noted that: 
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“Generally speaking, the classical limit of a quantum state is not a single 
classical orbit, but an ensemble of orbits. ” 

Moreover, Huang [17], when using Wigner distribution function [4], compared with a 
given probabilistic dehnition of classical pure ensemble, concluding that: 

“A wave function does not describe an individual particle but a classical pure 
ensemble. ” 

But not all of the classical limits prevent an individual classical particle description. 
For simplicity, we restrict our study to a one dimensional bound potential of U{x). As 
an exercise, we can show that in this restricted condition, there are no degenerate wave 
functions [18], so that every energy level En belongs to a wave function Now, we 

classify these wave functions based on pre and post measured states. 

1- Before-measurement wave function: A wave function before measurement 
is assumed to be represented by a superposition of all the possible eigenstates. 

CX) 

fj{x) =^Cn'llJn{x), (26) 

n =0 

Clearly, n —)■ cx) limit for such a wave function is meaningless. However, we can 
correspond a Pi{x) to every in terms of distinct energies {Eqm = E^). Roughly 

2 1 

(x)| —^ Pi{x) oc Pi{x) ^ Pj{x) Ei ^ Ej. 

^J[Ei = E) -U[x) 

Then, the set of {po{x),pi{x),p 2 {x),...} would describe an ensemble of non-identical 
particles. These particles would be non-interacting in case it is assumed that 
independent terms in a linear combination such (26) do not interact with each other. 

2- After-measurement wave function: On the other hand, when the 
measurement is performed, the system is thrown into one of the eigenstates [19] of 
Ri. In other words 

'lp{x) = lfn{x), 

then we can use n —)> oo limit for this 'if{x) to correspond it to a specihc p{x) (in a given 
energy Ed = En) and describe the probability density of an individual classical particle 
or an ensemble of identical particles. 

The difference between the classical limit of a wave function in a linear combination 
of independent eigenfunctions ||, and the classical limit of one of those independent 
eigenfunctions (here, a non-degenerate energy eigenstate) could be important when 
encountering various kinds of classical limits through literature. 

II Special examples for eager student are equation 18 in [6] and equation 19 in [17]; however the second 
one (free particle) has degeneracy. 
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6. Conclusion 

In this note, we narrowly focused on where a quantum feature of “dependence on 
the number of states at low energy levels” seemingly was lost when evaluating the 
dimensionless uncertainty of bouncing ball and harmonic oscillator potentials, since the 
quantum results were exactly the same as the classical results. We concluded that 
this loss was removable with changing potential intervals in terms of parity. In short, 
this meant including even wave functions in the solutions of the linear potential and 
excluding them for the quadratic potential by changing potential intervals. Hence it is 
concluded that the existence of even and odd parity solutions of Schrodinger equation 
is an inherently quantum mechanical property that is pursuable in the dimensionless 
averaging of lower energy levels, even in potentials for which Ehrenfest theorem is exact. 
This could be important, since dimensionless analysis must not eliminate inherently 
quantum features in general. 

Calculations also contain higher powers of position for potential to expand 
dimensionless analysis to other examples via the numerical method. According to figures 
1 and 2, higher powers of potentials depend on n in both symmetric and non-symmetric 
forms, while linear and quadratic potentials lose their n-dependence in the shape of 
bouncing ball and harmonic oscillator potentials, respectively. Therefore, we have not 
paid more attention to higher powers, where n-dependence can be explicitly seen in 
figures 1 and 2. Students could study our considerations for even larger powers of b, 
where they eventually might encounter numerical constrains. It would be an instructive 
challenge to hnd sources for lack of accuracy (whether or not correlated with execution 
time) in their computer codes to develop numerical skills. 

As other potential examples to study dimensionless uncertainty and compare it 
with classical dispersion, we can consider supersymmetric partner potentials for the 
potentials from the power of coordinate {U{x) oc x’^). A review of supersymmetry in 
quantum mechanics is explained in the article [20] or in the book [21]. In the case of 
quadratic potential, the supersymmetric partner potentials do not change, except for a 
constant. Therefore, wave functions fjnix), classical density p{x) and related Xmax do not 
change, so that averaging on dimensionless position remains unchanged, unlike averaging 
on dimensionless momentum that changes according to the changed Pmax, which causes 
separate levels in the amount of uncertainty (or even classical dispersion) in terms of 
quantum number n, just as separate levels of energy. Along with this special case, 
supplementary figure data is provided for the uncertainty of supersymmetric partner 
potentials from powers of coordinate, illustrating this separation and comparing it with 
figure 1. 

Finally, in the discussion section, we mentioned that approaching quantum 
averaging to that of classical in large quantum numbers does not vanish quantum 
uncertainty in the classical limit. Also, in the case of non-degenerate states, we have 
emphasized that it is the classical limit of a “before-measurement wave function” that 
can not describe an individual classical particle, while classical limit of an “after- 
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measurement wave function” (a non-degenerate energy eigenstate 'ifni^)) can describe 
an individual classical particle or an ensemble of identical particles. For a more 
complete study, we can look forward to extend the classification in terms of every linear 
combination of independent eigenfunctions (including degenerate states, rather than just 
a before-measurement wave function) versus one of those independent eigenfunctions. 

One pedagogical review of dimensional analysis has been explained in a recent 
paper [22] with a variety of examples. Pedagogically speaking, dimensionless analysis, 
as a secondary branch of dimensional analysis, can be given more attention because 
of its beneht for analytic and numerical problems. Two helpful examples of using 
dimensionless analysis are deletion of h in uncertainty and comparison with classical 
dispersion possible, and development of differential equations as dimensionless ones. 
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Appendix A. Derivation of phase space averages from integrating on 
position space 

The position probability function is obtained by integrating over the momentum variable 
P 


p{x) = / dpp(x,p) = constant. dp 6 




E 


(A.l) 


By using the properties 6{ax) = 5{x)/\a\ and 5{x'^ — of) = [(5(x -1- a) -|- 5{x — a)]/2|a| of 
the Dirac delta function, the classical probability distribution is reduced to 


p{x) = constant. / dp2m(5 (p^ -|- 2m[U(x) — Ej) = constant ..i 


2m 


|E - U(x)| 


X 


dp 


5 fp + 


+5 (p - V2m[D(x) - Ej) oc ^ (A.2) 

, which is in agreement with (2). In the next step, the phase space averages of any 
arbitrary function f{x,p) of position and momentum variables get reduced to those 
evaluated with the position probability distribution function p{x) as follows: 

if ix,p))= f dx f dpp{x,p)f{x,p) (X f dx f dp5 (^ + U{x) - e\ f{x,p) 
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oc J 

oc J dxi 

so that 


2m 


E - U{x) 

2m 
E - U{x) L 

1 


dp 5 + \/2m[U{x) — E]^ + ^P — ^/2m[U{x) — E]^ 


fix,-^2m[E-Uix)]) + fix, ^2m[E-Uix)]) 


(fix, P)) = 2 / -V2m[E- t/(x)]j + / ^/2m[E^^TJi^] 


fix,v) 

(A.3) 

(A.4) 


Appendix B. Matrix representation structure for solving Schrodinger 
equation 

In the numerical solutions, we are looking for eigenvalues and eigenvectors of the 
Hamiltonian operator 

d^ 


H = 


Uix), 


(B.l) 


2m dx'^ 

in its matrix representation. Suppose that Xj+i = Xi + Ax. For a given small Ax, we 
can write 


fU.) ~ /(^»+l) - fjXi) ^ _ fiXi+i)-fiXi) 

* Ax ’ * Ax 

fixi+ 2 ) - 2/(xi+i) + fixi) 


fix.) ^ 


Ax"^ 


(B.2) 

(B.3) 


In Dirac representation, | /) is a column matrix in the x basis 

/ fixi) \ 

fix2) 

I/) ^ \ 

\ fiXn) J 

For derivative operator D, we have 

ci/) = 1 /'). 


(B.4) 


(B.5) 


or 




( fixi) \ 

( 

D 


f{x2) 

1 



Ax 


nxn 

V fiXn) ) 

nx 1 ' 


/ fix 2 ) - fixi) \ 
fixs) - fix 2 ) 


(B.6) 


Therefore 


D = 


Ax 


( -1 1 
-1 1 0 
-1 1 

0 -1 1 
V 




(B.7) 


• ••• / 
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In a similar way, for the second derivative, we can write 

/ 1 -2 1 


= 


Ax 




1 -2 1 
1 -2 1 
0 1-2 


0 


V 


(B.8) 


/ 


According to the dehnition of a derivative, however, this matrix conld be a little different. 

Having written potential arrays in the shape of a diagonal matrix, we can have 
a matrix representation of Hamiltonian. Generally speaking, making a dimensionless 
Schrodinger equation is more formal at first, but for our purpose, which is averaging on 
dimensionless quantities of (8), multipliers like m and h in the Hamiltonian and also, a 
in the potentials (14) and (15) can be ignored easily. 

Note that wave functions naturally tend to zero at x —)> oo, but since writing 
a matrix representation with inhnite arrays is impossible, then our wave functions 
(eigenvectors) do not continue to infinity, so we have to adopt an end for x axis like A. 

X G [—A, A] for symmetric potentials, 

X G [0, A] for non-symmetric potentials. (B.9) 

For a bound state, A can be chosen about a few times more than classical return points 
by considering enough Xj in the interval. 

Appendix C. Moments for half wave functions of symmetric harmonic 
oscillator 


Although odd wave functions of the symmetric harmonic oscillator potential are 
solutions of non-symmetric harmonic oscillator (21), we evaluate moments for either 
even or odd wave functions by starting from zero up to inhnity, which we call half wave 
functions, for pedagogical goals or possible implications. The half wave functions of 
symmetric harmonic oscillator are 


\i^n{x)f = + gxp 

TT22'^~^n\ 


-(2n + l) 


X 




(2n + 1)2 


1 X 


A. 


(C.l) 


, where An satishes (n-|-l/2)hci; = l/2mci;^A„^. Note that wave functions are normalized 
in the interval [0, cxd). The moments are as follows 


k 1 

(^^)n = f [Hn{y)?dy, y = X 


or 


7r22"- ^n\ 




Ar. 


(2n + l)-| 
7r^2"“^n! 


1 /^exp (-1/2) [Hn{y)fdy. 


(C.2) 


(C.3) 
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There is a useful tabled integral as [23] 

OO 

r [(m + l)/2] 


exp(—Q!a;^)da: = 


r) 

2a 2 


16 


(C.4) 


, where r(x) indicates the Gamma function. The moments can be written by a 
summation of these integrals. In this way, mean values can be obtained. However, 
it may not simply give us a general term in respect to n for (C.3). Instead, it can be 
seen that for even numbers k, integrand is an even function. Hence 




(2n+ 1)- 
7r^2^n! 


2/^™exp(-|/^) [Hn (?/)] dy, k = 2m, meN 


(C.5) 


They can be evaluated using orthogonality and recurrence relations for Hermite 
functions (see, for example, exercise 13.1.11 of [24]) and substituting recurrence relation 
again for even levels higher than m = 1. For instance. 


/ 4, 3n2 + n+l/2 

/n - T 


lim (X^) = 

\ / n Q 


2 (2n + 1)^ ’ 

Here we return to the first moment by using derivative formulas [23] 

^ [exp{-x‘^)Hn (x)] = - exp{-x‘^)Hn+i (x ), (x) 

and integral recurrence relation 

OO OO 

/ xexp (—[Hn (x)]^ dx = - exp{—x^)Hn+i (x) Hn {x) dx 


(C.6) 


(C.7) 


+n / ex'^{—x^)Hn{x) Hn-i{x) dx. 


(C.8) 


Applying derivative formulas and then integrating by parts, we have 

OO ^ 

J xexpi-x"^) [Hn{x)f dx = ^ [Hn{0)f + n\^ {n^-k)\ 

n k=l 


1 "" 






Remembering 


H2n (0) = (-!)' 


(2n)! 


n\ 


, H2n+1 (0) — 0 


and 


E 


[2 {n-k)]\^2k ^ ^ri + 1)! 


[{n-k)\y 


n! 


2 ’ 


(C.9) 


(C.IO) 


(C.ll) 
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we can derive 


OO 

J xexp{-x‘^)[H2n (x)]^da; 
0 


1 

2 



+ (2n)!5^ 

k=0 


[2 (^- fc)]! o2fc 

[{n-k)\f 


1 

2 


■(2n)!l^ 

n\ 


+ (2n + l) 


'{2n)\V 

n\ 



[H2n (0)] = 


(C.12) 


and 


a:exp(-a:^) [ifsn+i (a:)]^da; = (2n + 1)! ^ 

k=o 


2(2n + l)!5^ 


[2 


^0 


= 2 


(2n + 1)! 


k=0 
2 


n\ 


= 2{2n + lf[H2n ( 0 )]' 


So, the first moment leads to 

2 (2n + l) (2n)! 


{X)2n 

{^)2n+l 


7r2 (477, + 1)2 22"- {n\f 

2 (2n + l) (2n)! 


vr2 (477 + 3)2 22^ (tt!) 

The limit of n —)■ cx) can be evaluated by using ^ 

(2?7,)! 1 

22 -( n !)2 ^ 7 ^’ 


(C.13) 


(C.14) 

(C.15) 


(C.16) 


thus 


lim {X) 

n—^oo 


2 n 


lim {X) 

n—>-oo 


2n+l 


2 77.2 (2?7,)! 

—j- lim -^ 

yj-2 n-7oo 22”' (77!) 


2 

71 


(C.17) 


Equations (C.6) and (C.17) can be checked against the classical results when we average 
over a quadratic classical density function. 


^ Equations (C.ll) and (C.16) have been obtained by MAPLE mathematical software. 
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Supplementary Figure 



b 

Caption Cl. Classical (dashed) and quantum (solid) uncertainty of dimensionless 
variables for supersymmetric partner potentials of equation (14)(6 = 1, 2, 3,4) on three 
first levels 









